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Abstract 

We use Gelfand-Tsetlin diagrams to write down the weight multiphcity function for the Lie algebra 
slfeC (type Ak-i) as a single partition function. This allows us to apply known results about partition 
functions to derive interesting properties of the weight diagrams. We relate this description to that of 
the Duistermaat-Heckman measure from symplectic geometry, which gives a large-scale limit way to 
look at multiplicity diagrams. We also provide an explanation for why the weight polynomials in the 
boundary regions of the weight diagrams exhibit a number of linear factors. Using symplectic geometry, 
we prove that the partition of the permutahedron into domains of polynomiality of the Duistermaat- 
Heckman function is the same as that for the weight multiplicity function, and give an elementary proof 
of this forsUC (^3). 



1 Introduction 

For a long time there has been a lot of interest, both in mathematics and physics, in finding ways to 
determine with what multiplicity a weight appears in the weight space decomposition of a finite-dimensional 
irreducible representation of a semisimple Lie algebra. Although there is a multitude of formulas to compute 
these multiplicities, involving partition functions (Kostant's formula), recursions (Freudenthal's formula), 
counting paths (Littelmann's formula), this is still a computationally hard problem. For type A (SL^C, 
GLfcC, SU(fc)), these multiplicities are known to be the Kostka numbers, which express the Schur symmetric 
functions in terms of the monomial symmetric functions. 

Here we explore the structure of the weight diagrams in type A, not from a symmetric functions perspec- 
tive, but using an array of tools from combinatorics, convex geometry and symplectic geometry, such as 
Gelfand-Tsetlin diagrams, Kostant's multiplicity formula, and the so-called "Quantization Commutes with 
Reduction" Theorem. We describe how the weight diagrams are partitioned into domains of polynomiality, 
and how this is related to the Duistermaat-Heckman function studied by symplectic geometers. 

After a brief reminder about the structure of the Lie algebra sl^C, we introduce our main tools, Gelfand- 
Tsetlin diagrams and partition functions. Gelfand-Tsetlin theory provides a way of computing weight 
multiplicities by counting certain combinatorial diagrams, or equivalently, by counting the number of in- 
teger lattice points inside certain polytopes. We will use this and some notions from linear and integer 
programming to reduce this counting problem to evaluating a single partition function. 

Theorem EH For every k, we can find integer matrices Ek and Bk such that the multiplicity function for 
sIfcC can he written as 

' X 



mx{l3) = (pEk Bk 



Expressing the multiplicities as a single partition function allows us to use general facts about partition 
functions and their chamber complexes to derive interesting properties of the weight diagrams. For example, 
the multiplicities have the following polynomiality property. 
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1 Introduction 



Theorem l5.1l There is a chamber complex C^'^^ on which the weight multiplicity function is determined by 
polynomials of degree ('^2^) in the Pi, with coefficients of degree C^j^) in the Xj. 

From this theorem we can deduce a pointwise scaUng property (i.e. for fixed A and P). This property 
fCorollarv l5.2|l was known already in the context of symmetric function theory, where it was proved using 
a fermionic formula for the Kostka-Foulkes polynomials (see HOI). It shows that although the Gelfand- 
Tsetlin polytopes are not always integral polytopes |24| , their Ehrhart quasipolynomials are in fact always 
polynomials. 

The partition of the weight diagram into its domains of polynomiality can be described explicitly. The 
convex hull of a weight diagram is a permutahedron. There is in symplectic geometry a function on the 
permutahedron, called the Duistermaat-Heckman function, that approximates the weight multiplicities 
and is known to be piecewise polynomial. Its domains of polynomiality are convex subpolytopes of the 
permutahedron, and there is an explicit description of the partition in terms of walls separating the do- 
mains. Using known results on quantization and reduction of symplectic manifolds, we can prove that the 
Duistermaat-Heckman function and the weight multiplicity function give rise to the same partition of the 
permutahedron. 

Theorem 13.21 The partitions of the permutahedron for su(k) (or slkC) into its domains of polynomiality 
for the weight multiplicities and for the Duistermaat-Heckman measure are the same. Namely, the walls are 
determined by convex hulls of the form conv(W • cr(A)) where a G 6k and W is any parabolic subgroup of 
&k generated by all reflections corresponding to roots orthogonal to a conjugate of a fundamental weight. 

In Kostant's multiplicity formula, multiplicities are expressed as a sum of partition functions evaluated at 
k\ points shifted by a factor depending on the choice of a positive root system. We can take advantage of 
the apparent lack of symmetry of Kostant's multiplicity formula to find interesting factorization patterns 
in the weight polynomials of the boundary regions of the weight diagrams. 

Theorem 16.21 Let R be a domain of polynomiality for the weight diagram of the irreducible representation 
of slkC with highest weight X, and pr be its weight polynomial. Suppose that R has a facet lying on the 
boundary of the permutahedron for X that has Oiujj) as its normal vector, for some 9 G &k. If "f = 7(A) is 
the defining equation of the hyperplane supporting that facet, then pn is divisible by the j{k — j) — 1 linear 
factors 7 + 1, 7 + 2, . . 7 -|- j(/c - or 7 - 1, 7 - 2, . . "f-jik-j) + 1. 

The main tool for proving this theorem is a family of hyperplane arrangements, called Kostant arrangements, 
on whose regions we have different polynomials giving the multiplicities. The Kostant arrangement also 
provides a method for finding linear factors in the difference between the weight polynomials of two adjacent 
regions. A generalization of the Kostant arrangements is also essential to the proof of Theorem 15. II which 
establishes that although in general we get quasipolynomials in the chambers of the complex associated to 
a vector partition function, we get polynomials for the weight multiplicity function in type A. 

Theorem 16.51 Let Ri and R2 be two adjacent top- dimensional domains of polynomiality of the permuta- 
hedron for a generic dominant weight X of slkC, and suppose that the normal to their touching facets is 
in the direction cr^LUj) for some a G 6fe. If pi and p2 are the weight polynomials of Ri and R2, and 7 is 
the linear functional defining the wall separating them, then the jump pi — p2 either vanishes or has the 
j{k — j) — 1 linear factors 

(7 - + 1), (7 - + 2), . . . , 7, ■ • ■ , (7 + s+ - 2), (7 + s+ - 1) 
for some integers s^ , s^ > satisfying 

s^ +s+ ^ j{k~j). 

Similar factorization phenomena were recently observed to hold for general vector partition functions by 
Szenes and Vergne |5T] . 

Finally, we explicitly compute the chamber complex for ^3, and find it is not optimal, but that we can glue 
together parts of it to obtain a simpler complex. We can deduce symbolically from the form of this complex 
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that the optimal partitions of the permutahedron for A3 under the weights and the Duistermaat-Heckman 
measure are the same. Computing the chamber complex for A3 is nontrivial because of the complexity 
of the arrangement. To the best of our knowledge, these computations for generic dominant weights of 
A3 have not been done. A study was done by Guillemin, Lerman and Sternberg in |14| for some of the 
degenerate cases when A has a nontrivial stabilizer. The number of domains of polynomiality turns out to 
be significantly larger than they originally suspected. 

1.1 The Lie algebra sl^C (type A^^i) 

The simple Lie algebra sf^C is the subalgebra of gl^C ~ End(C'^) consisting of traceless k x k matrices 
over C. We will take as its Cartan subalgebra f) its subspace of traceless diagonal matrices. The roots and 
weights live in the dual f)* of f), which can be identified with the subspace + • • ■ + = of R*^. The roots 
are {ci — Cj : 1 < i j < k], and we will choose the positive ones to be A+ = {ci — ej : 1 < i < j < k}. 
The simple roots are then = — ej+i, for 1 < i < k — 1, and for these simple roots, the fundamental 
weights are 

oji = — (fc — i, k — i, . . . , k — i, —i, —i, . . . , —i) , 1 < i < ^ ^ 1 • (1) 
k " . ' « ' 

i times k — i times 

The fundamental weights are defined such that {ai.uuj) = 5ij, where (■, •) is the usual dot product. The 
integral span of the simple roots and the fundamental weights are the root lattice A/j and the weight lattice 
h-w respectively. The root lattice is a finite index sublattice of the weight lattice, with index k — 1. 

For our choice of positive roots, 

k—1 

E "-E^fc = ^(^-l'^-3,...,-(fc-3),-(fc~l)). (2) 

The Weyl group for s[feC is the symmetric group Sfe acting on {ei, . . . , e^} (i.e. cr(ei) = &a(i))i f^nd with the 
choice of positive roots we made, the fundamental Weyl chamber will be Co — {(Ai, . . . , A^) : ~ 
and Ai > • • • > Afc}. The action of the Weyl group preserves the root and weight lattices. The Weyl 
orbit of a weight A is the set Sfc • A — {cr(A) : a £ 6k}- We refer to the convex hull of &k ■ A as the 
permutahedron associated to A. Weights lying in the fundamental Weyl chamber are called dominant, and 
we will call elements of the Weyl orbits of the fundamentals weights conjugates of fundamental weights. 

The finite dimensional representations of slfeC are indexed by the dominant weights Aw n Cq, and for a 
given dominant weight A, there is a unique irreducible representation p\ : sI^C — > ^[(Va) with highest 
weight A, up to isomorphism. Details about their construction are well-known and can be found in |1(J | or 
|11| . for example. We have the weight space decomposition according to the action of f) 

^a = 0(^a),. (3) 

The weights of this representation (those /3's for which (Va)^ ^ 0) are finite in number, and they can be 

characterized as follows (see JHl): they are exactly the points /3 of the weight lattice Aw that lie within 
the convex hull of the orbit of A under the Weyl group action, denoted conv((3fc • A), and such that X — (3 
lies in the root lattice. Hence 

V,^ (^a),. (4) 

/3e(A+AH)nconv(efc-A) 

The multiplicity m\{[3) of the weight (i in V\ is the dimension of (Va)^7 and all the conjugates of (3 under 
&k have the same multiplicity. The weight diagram of V\ consists of the weights of V\ (as a subset of Aw) 
together with the data of their multiplicities. 
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1.2 Gelfand-Tsetlin diagrams 



There are several ways to compute weight multipHcities. An important one is Kostant's multipUcity formula 
|21j . which can be deduced from Weyl's character formula (see We first need to define the Kostant 

partition function given a choice of positive root system A_|_ : 



K 



(5) 



i.e. K(v) is the number of ways that v € [}* can be written as a sum of positive roots. 
Kostant's multiplicity formula [2] is then 

mxW - i-^y^^^KiaiX + 6) -{(3 + S)) , 



(6) 



where l{a) is the number of inversions a. Kostant's partition function and multiplicity formula extend to 
all complex semisimple Lie algebras. See jlSi for more details. 



1.2 Gelfand-Tsetlin diagrams 



Gelfand-Tsetlin diagrams were introduced by Gelfand and Tsetlin jl2j as a way to index the one-dimensional 
subspaces of the (polynomial) representations of GL^C Their construction relies on a theorem of Weyl 
that describes how the restriction to GLj;_iC of an irreducible representation of GL/jC breaks down into 
irreducible representations of GL^.iC (see jniEll^)- They are equivalent to semistandard tableaux (see 
|13p. but they have a "linear" structure that we will exploit. 

Definition 1.1 Let v — [vi, . . . , Vm) o,nd 7 — (71, . . . , 7m-i) be two partitions. We will say that 7 interlaces 
ly, and write j <iv, if 

1^1 > 71 > 1^2 > 72 > 1^3 > ■ ■ ■ > Vm-l > 7m-l > ■ 

Theorem 1.2 (Weyl's branciiing rule | 13L VA'6\ ) Let p\ he the (polynomial) irreducible representation 
0/ GLfeC with highest weight A = Ai > A2 > • . . A^ > 0. The decomposition of the restriction of p\ to 
GLfe_iC into irreducible representations o/GLfe_iC is given by 

PAi =0Pp- (7) 

After restricting p\ to GLfc_iC and breaking it into GLfc_iC-irreducibles, we can restrict to GLfc_2C: 



PA I 




Again, we can apply Weyl's branching rule to each to break them into irreducible representations of 
GLfe_2C to get 

PAi = P.- (9) 



We can keep going recursively, and for convenience, let us denote by A^™^ — A^™^ > • ■ ■ > aIT^ > the 
partitions indexing the irreducible representations of GL.,„C. We then get that 



PA I 



'<=LiC A(i)<l---<A('')=A 



Pa(i) ■ (10) 
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Definition 1.3 A sequence of partitions of the form A^^-' < • • ■ <l A^^*^') ~ \ is called a Gelfand-Tsetlin 
diagram for A, and can he viewed schematically as 



(11) 



,(fe) ,(fc) ,(fc) 

^1 '^2 -^k-l 

,(fe-l) ,(fe-l) 

Af) A^^) 

with \^^^ — Xj and each Aj*"* is a nonnegative integer satisfying 

^+1) > > ^(i+i) for l<i<k-l,l<] <i. (12) 

Let be the one-dimensional subspace of V\ corresponding to a Gelfand-Tsetlin diagram T). It is shown in 
1331 that V-D lies completely within one weight space in the weight space decomposition of V\: Vd C (Va)^ 
if 

m rn— 1 

/?„ = ^ aI'") - A^"-'^ for 1 < TO < fc (13) 

i=l i=l 

or, equivalently, 

m 

/3i + • • ■ + /3m = 51 for 1 < m < fc. (14) 

Hence Gelfand-Tsetlin diagrams for A correspond to the same weight if all their row sums are the same. 
This discussion is summarized in the following theorem due to Gelfand, Tsetlin and Zelobenko. 

Theorem 1.4 f | 12L I33 p For A = (Ai, . . . , Afc), the number of Gelfand-Tsetlin diagrams with first row A 
is the dimension of the irreducible representation V\ of GL^C with highest weight A. Furthermore, the 
multiplicity mx{/3) of the weight (3 in the irreducible representation of Ghk'C with highest weight A is given 
by the number of Gelfand-Tsetlin diagrams with first row A such that equation (or ) is satisfied. 

Two irreducible representations V\ and of 0[j,C restrict to the same irreducible representation of slkC 
if Ai — 7i is some constant independent of i for all i. Hence we might as well require that the A^ sum 
up to zero. However, normalizing the sum this way can introduce fractional values of A, so we'll have to 
translate A back to integer values when writing down Gelfand-Tsetlin diagrams for those representations, 
or, equivalently, translate the integer lattice along with A, so that the inequalities 

> ^ > fQj. l<i<k-l,l<j <i, 

always have 

a5.'+'^ - A**) e N and A^ - \^+^ e N . 

There is a geometrical way to view the enumeration of the number of Gelfand-Tsetlin diagrams for a given 
A. With AC^) = A fixed, we can let all the other variables {A^™'' : 1 < i < to, 1 < to < fc} be real variables. 
The system of inequalities (|12f) among the entries of Gelfand-Tsetlin diagrams, when viewed over the reals, 
defines a rational polytope, called the Gelfand-Tsetlin polytope for A and denoted GT^. GT^ has dimension 
at most (2), and equal to that number if the A^'s are distinct. We can consider the intersection of this 
polytope with the affine subspace obtained by fixing a weight f3 (fixing the row sums using equations H13(l 
or (|14|l ). We also get a rational polytope this way, called the Gelfand-Tsetlin polytope for A and [3 and 
denoted GTa./j . Its dimension is at most ('' j^)- Kirillov conjectured in that the polytopes GT^,/? are 
integral polytopes, but this was recently disproved by De Loera and McAllister 24 . 

The upshot is that integer solutions to the Gelfand-Tsetlin diagram constraints then translate into integer 
points inside the polytopes, hence the number of Gelfand-Tsetlin diagrams of weight (3 for A is the number 
of integer points in the polytope GT^^^. 
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1.3 Partition functions and chamber complexes 

Partition functions arise in the representation theory of the semisimple Lie algebras through Kostant's 
formula for the multiplicities Kostant's partition function sends a vector in the root lattice to the 
number of ways it can be written down as a linear combination with nonnegative integer coefRcients of 
the positive roots, and this is a simple example of a more general class of functions, called vector partition 
functions. 

Definition 1.5 Let M be a dxn matrix over the integers, such that kerMnM"Q — 0. The vector partition 
function (or simply partition functionj associated to M is the function 

0M : Z'^ — > N 

b ^ |{xeN" : Mx^b}\ 

The condition kerA/ n R"q — forces the set {x £ N" : Mx — b} to have finite size, or equivalently, the 
set {x £ M"q : Mx = 6} to be compact, in which case it is a polytope Pb, and the partition function is 
the number of integral points (lattice points) inside it. 

Also, if we let Mi, . . . , Mn denote the columns of AI (as column- vectors), and x = {xi, . . . , Xn) G ^>oi then 
Mx = xiMi + 2:2-^2 + • • • + XnMn and for this to be equal to b, b has to lie in the cone pos(M) spanned 
by the vectors Mi. So ipM vanishes outside of pos(Af). 

It is well-known that partition functions are piecewise quasipolynomial, and that the domains of quasipoly- 
nomiality form a complex of convex polyhedral cones, called the chamber complex. Sturmfels gives a very 
clear explanation in [30] of this phenomenon. The explicit description of the chamber complex is due to 
Alekseevskaya, Gel'fand and Zelevinskii There is a special class of matrices for which partition functions 
take a much simpler form. Call an integer dxn matrix M of full rank d unimodular if every nonsingular 
d X d submatrix has determinant ±1. For unimodular matrices, the chamber complex determines domains 
of polynomiality instead of quasipolynomiality 30 . 

It is useful for what follows to describe how to obtain the chamber complex of a partition function. Let 
M be a d X n integer matrix of full rank d and (pM its associated partition function. For any subset 
a C {1, . . . denote by M^- the submatrix of M with column set ct, and let = pos{M^), the cone 
spanned by the columns of Mo- . Define the set B of bases of M to be 

B = {a C {1, . . . ,n} : |(t| = d and rank(Mo) = d} . 

B indexes the invertible dx d submatrices of M. The chamber complex of (j)M is the common refinement of 
all the cones t^, as ct ranges over B (see ^). A theorem of Sturmfels [301 describes exactly how partition 
functions are quasipolynomial over the chambers of that complex. 

1.4 The chamber complex for the Kostant partition function 

If we let Ma„ be the matrix whose columns are the positive roots A^"-* of An, written in the basis of 
simple roots, then we can write Kostant's partition function in the matrix form defined above as 

Consider for example the simple Lie algebra sUC, or A3. The positive roots are A^"^'' — {ei — Cj : I < i < 

j < 4}. Writing the positive roots in the basis of simple roots, we have A^'''' — {ai, 02, as, ai + a2, 012 + 
as, ai + a2 + a^}. This gives 

/ 1 1 1 \ 
M^3 =010111 

\ 1 1 1 / 



1.4 The chamber complex for the Kostant partition function 
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which has the bases 

B = {123, 125, 126, 134, 135, 136, 145, 146, 234, 236, 245, 246, 256, 345, 356, 456} , 
where we're writmg 1112*3 for {ji,i27*3}- 

All the cones corresponding to these bases are contained in the first cone with basis {1, 2, 3} which is just 
the positive octant in R^. To picture the chamber complex, we can look at the intersection of these cones 
with the hyperplane x + y + z = 1. Figure^shows the cones given by the bases of B, while Figure |21 shows 
their common refinement (this originally appeared in j25) ) . Finally, since it is readily checked that Ma^ 
is unimodular, this shows that the Kostant partition function for A3 has 7 domains of polynomiality. It 
is an open problem mentioned by Kirillov in |2U| to determine the numbers of chambers for the Kostant 
partition functions for the Lie algebras A^- De Loera and Sturmfels 2^ have computed the numbers for 
n < 6 and computed the polynomial associated to each chamber for n < 5. 




123 125 126 134 135 136 145 146 




234 236 245 246 256 345 356 456 



Figure 1: Basis cones for the Kostant partition function of A3. 



"2 




Oil 013 



Figure 2: Chamber complex for the Kostant partition function of ^.3. 

The following lemma is a well-known fact about Ma„ and can be deduced from general results on matrices 
with columns of O's and I's where the I's come in a consecutive block (see \27\). 

Lemma 1.6 The matrix Ma„ is unimodular for all n. 

Ma„ unimodular means that the Kostant partition functions for A„ is polynomial instead of quasipolyno- 
mial on the cells of the chamber complex. In general, for AI unimodular, the polynomial pieces have degree 
at most the number of columns of the matrix minus its rank (see ^30,). In our case, Ma„ has rank n and 
as many columns as An has positive roots, ("^^) • Hence the Kostant partition function for An is piecewise 
polynomial of degree at most ("J^) — n = (2). 

Remark 1.7 In view of Kostant 's formula for the weight multiplicities this means that the multiplicity 
function m\{j3) for An is piecewise polynomial of degree at most (2) in the [3 -coordinates if the \-coordinates 
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2 The multiplicity function as a single partition function 



are fixed, or degree (2) in the X-coordinates if the f3- coordinates are fixed. So we can regard it as a piecewise 
polynomial function of degree Q) in the j3i 's, with coefficients of degree (2) in the Xj 's. This will be made 
precise in Sections^ and\^ 

From now on, we will be interested in the multiplicity function for slfeC, of type Ak-i, and thus use the 
results above with n — k ~ I. 



2 The multiplicity function as a single partition function 

Our first theorem presents a new conceptual approach to computing multiplicities. This approach is efficient 
for large A in low ranks, ft has the additional advantages of allowing us to use known facts about partitions 
functions. 



Theorem 2.1 For every k, we can find integer matrices and _Bj. such that the multiplicity function for 
slkC can be written as 

"^a(/3) = 0£. (fifef ^ )V (15) 



Proof. Consider a Gelfand-Tsetlin diagram where we will think of A = (Ai, . . . , Afe) and f3 = . . . , Pk) 
as parameters, with the conditions that — ft — 0- The variables in the diagram are A^*'' with 

l<i<fc— 1, I < j < i- Each of these (2) variables is wedged between two entries at the level above, so 
we get a system of 2(2) = k{k — 1) inequalities. Using equation p4l) . relating the row sums to the /3i's, we 



can get rid of the k — 1 variables A^^'', A2^\ 



A 



(fc-i) 
fc-i 



Ai 



A^ 



(fc-i) 



>(fc-i) 



Afe-1 



(/3i 



-Pk = 0) 
-fik-i) 



(3) 



i(2) 



,(3) 



.(2) 



v(3) 



+ P2 + 133) 



(16) 



The remaining variables (boxed in the above diagram) are A^*^ with l<i<k — 1, l<j<i — 1 and there 
are C^j^) them. To get a system in partition function form, we need to transform the inequalities into 
equalities satisfied by nonnegative variables, however the A^*^ can take negative values. Let 



l<i<fc-l, l<i<i-l 



be the differences between the variables and the ones immediately above and to the right of them, recalling 
that A^'^'' = Xj. Upon doing the substitution in the system of inequalities, {'^2^) of the inequalities simply 
become sj'' > because of equation (|12|) . So we are left with a system of = k{k — 1) — (''2^^) = 

^(fc — l)(fc + 2) inequalities in the K = {^^^) nonnegative and integral variables s^j \ which we will relabel 
si, . . . , sk for convenience. 



The final step is to transform the inequalities into equalities. To this effect, we write each inequality in the 
form 

k k 



amKSK < ^ bmjXj + ^ CmjPj , 
J=l 3 = 1 



2 The multiplicity function as a single partition function 
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for 1 < m < and integers Qmi, ■ • • , flmif, &mj, c^j (1 < J < fc)- 

We introduce a slack variable for each inequality to turn it into an equality: 

k k 

flrnlSl + am2S2 + ' ' ' + a„iKSK + SK+m = ^ &mj + ^ Cmjl3j ■ 

The slack variables s^+i, ■ ■ • , sk+n are nonnegative, just like the previous K Si, and integral solutions to the 
system of inequalities will correspond to integral solutions to this system of equalities, so sk+i, . . . , sk+n 
are not only nonnegative but integral. 



Finally, we can write the system of equalities in matrix form: 



/an • • • aiK 
\ ajvi • ■ • aNK 



IN 



Eh 



SK 
SK+1 

\ Sk+n J 



(17) 



The result follows, since the number m\{(3) of integral solutions to the Gelfand-Tsetlin inequalities is the 
number of all integral nonnegative solutions to this matrix system. ■ 

The partition function (pEk in the above theorem lives on a larger dimensional space than the one we need. 
It takes values in R-^ = M('=~i)('=+2)/2^ whereas the part that interests us, the space given by as the 

\i and (3j range over M, has dimension 2k — 2. Let 



B 



Bk 



k k \ 

i=l i=l ) 



(18) 



then the only part of the chamber complex that is relevant to the multiplicity function is its intersection 
with B. Since the chamber complex is obtained as the common refinement of the base cones, we will 
get the same thing if we find the refinement of the base cones and then intersect the result with B, or 
intersect the base cones with B first and then find the common refinement of those restricted base cones. 
Since we only need the restricted chamber complex, this simplifies the computation because we have to 
deal with 2k — 2-dimensional cones instead of {k — l)(fc + 2)/2-dimensional ones. Another bonus we get 
from working on B is that on this space, Bk is an invertible transformation, so we can rectify the cones to 
(A, /3)-coordinatcs. In effect, we remove the coordinate "twist" due to matrix Bk- 



Definition 2.2 We will denote by C^'^^ this rectified {2k — 2) -dimensional complex in [\^ j3)- coordinates. 



Because Ek is not unimodular in general, the associated partition function will be quasipolynomial on the 
cells of the chamber complex. We will prove in Section that it is actually polynomial on the cells of 
the complex. As such, we will from now on refer to the domains of quasipolynomiality of the multiplicity 
function as domains of polynomiality. 

Remark 2.3 The multiplicity function also satisfies another sort of polynomiality property. There are 
many ways to think of fixed type A dominant weights A and (3 as living in st^C for any sufficiently large r. 
It is known (see for example ^IJ,^) that if m^^\(3) is the multiplicity of (3 in the irreducible representation 
Vx ofslrC, then m'^^\f3) is given by a polynomial function in r, for r large enough. Bounds on the degree 



10 



3 Domains of polynomiality via Duistermaat-Heckman theory 



of this polynomial are also given. This result is shown to extend to the other classical groups 'Iff) and also 
the classical ajjine Kac-Moody algebras j^. In our investigation of the weight multiplicities, we instead fix 
the rank of the Lie algebra and study the polynomial dependence in the X and (3 variables. 

Definition 2.4 For every A in the fundamental Weyl chamber, let 

L(A) = {(Ai,...,Afc,/3i,...,/?fc) : ft G M} . (19) 
Note that this space is really (k — 1)- dimensional since (ij ~ 0. Define also the projection 

pa: (Ai,...,Afc,/3i,...,/3fe) ^ (Ai,...,Afe). (20) 

Remark 2.5 The intersection of C'^'^^ with L{\) will give domains of polynomiality for the weight diagram 
of the irreducible representation o/sIfeC with highest weight A. The partition into domains that we get 
this way, however, is not optimal, as shown for sl^C in Section^ Some adjacent regions have the same 
weight polynomial and their union is again a convex polytope, so they can be glued together to yield a larger 
domain. 

Corollary 2.6 Let Cj^^ be the chamber complex given by the common refinement of the projections Pa(t) 
of the cones of C^'^^ onto M'^'. Then C^^^ classifies the X's, in the sense that if \ and A' belong to the same 
cell of Cjl^\ then all their domains are indexed by the same subsets of cones from C*-*^-*, and therefore have 
the same corresponding polynomials. 

Proof. We can index the top-dimensional domains by the top-dimensional cones t of CC'^). The domain 
indexed by cone r is present in the weight diagram (permutahedron) for A if and only if A £ pa{t). M 



3 Domains of polynomiality via Duistermaat-Heckman theory 

The chamber complex for the multiplicity function can be used to identify domains of polynomiality. 
However, these domains are not guaranteed to be as large as possible as seen in the examples of Section[71 In 
this section we improve the partition of the permutahedron into domains of polynomiality by identifying it as 
a bounded plane arrangement that appears in symplectic geometry. We begin by introducing the symplectic 
setup corresponding to the special case of type A multiplicities. Then we define the Duistermaat-Heckman 
function via an integral. This function is piecewise polynomial with natural domains of polynomiality in 
terms of Weyl group orbits. Finally we will use a powerful theorem of Meinrenken [221 and Vergne |32| . 
the so-called Quantization Commutes with Reduction Theorem, to show that the multiplicity function can 
be written locally as a very similar integral with the same domains. 

Let G = SU(fc), T the Cartan subgroup of G, q and t their Lie algebras, t;^ the fundamental Weyl chamber 
and Avi^ C t* the weight lattice of G. For A G fl Aw, we will denote by A a the convex hull of the Weyl 
group orbit of A in t* (i.e. the permutahedron associated to A). Let 0\~G- diag(A) be the coadjoint orbit 
for A. We can view 0\ as the set of fc x fc Hermitian matrices with eigenvalues {Ai, . . . , Afe}. By a theorem 
of Schur and Horn [281 117| (or Kostant's convexity theorem [33], which extends the result to all compact 
Lie groups), is the image of the coadjoint orbit Ox with respect to the projection map 

TT : 0* ^ t* . (21) 

The coadjoint orbits 0\ are the geometric counterpart to the irreducible representations of G with highest 
weight A. Note, the multiplicities for irreducible representations for SU{k) and SL{k) are the same. 

Consider M ~ 0\ and let $ : Af ^ t* be the restriction of tt to AI. In this case, $ is the moment map of 
the symplectic manifold M under the T action. The set Aiog C Aa of regular values of decomposes into 
a disjoint union of its connected components: 



Areg = U Ai 



(22) 
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and each is an open convex polytope by a generalization of Kostant's convexity theorem due to Atiyah 
P] and Guillemin-Sternberg ^^l- In fact, the singular values of <1> have the following nice combinatorial 
description. This theorem first appeared in Heckman's thesis ^B] 

Theorem 3.1 (\ 14l Theorem 5.2.1], [16 ) The singular points of the moment map $ : M ^ t* are the 

convex polytopes 

conv(VF • ct(A)) (23) 

where a G 6^ and W is any parabolic subgroup of &k generated by all reflections corresponding to roots 
orthogonal to a conjugate of a fundamental weight. 

In other words, the A^'s are the regions in the arrangement given by slicing the permutahedron by bounded 
hyperplane regions parallel to one of its exterior facets which pass through orbit points. See for example 
Figure 21 in Sectional Note, this is not a hyperplane arrangement inside a polytope since the convex hulls 
do not necessarily extend to the boundary of the permutahedron. 

Duistermaat and Heckman have shown that much of the geometry of coadjoint orbits can be determined 
simply by studying the Aj's. For e A^, the symplectic reduction of M at the regular value of $ is 
defined by 

M^ = <i>-\ti)/T. (24) 

For arbitrary G, the reduced space of M at a regular value of <i> is an orbifold, but for SU(fc) this 
orbifold is a compact Kahler manifold whose symplectic form we will denote hy 0. Duistermaat and 
Heckman [J] have shown that Af^ ^ -^^mo complex manifolds for any pair fio, ^J■ € A;. Furthermore, they 
have also shown the linear variation formula 0, 

= ^^Mo + (Ai - /-io,c) , (25) 

where c € 1(E) ri^(M^) is the Chern form of the principal T-bundlc <i>^^(/i) M^. Therefore, they use this 
fact about the symplectic forms to show that, for of dimension 2d, the symplectic volume function 

/fH(/i)=/ expc.^=/ (26) 

is a polynomial function on A.^ , called the Duistermaan-Heckman polynomial. Note that the only aspect of 
this integral that depends specifically on ^, and not just on which connected component of regular values 
contains it, is the symplectic form which is determined by 125|) . From the integral, one can show that the 
degrees of these polynomials are less than or equal to (dimM)/2 — dimG. 

Using a theory of quantization initiated by Kostant, Kirillov and Souriau (see |22j, for instance), we can 
apply the same reasoning used by Duistermaat and Heckman to the multiplicity function. 

Theorem 3.2 The partitions of the permutahedron for su{k) (or slkC) into its domains of polynomiality 
for the weight multiplicities and for the Duistermaat-Heckman measure are the same. Namely, the domains 
are the connected components of regular points determined by (|23|) . 

Proof. Let Td{M^g) be the Todd form of M^g. The Quantization Commutes with Reduction Theorem 
|26[ I32| asserts that for /i G A,;, 

mxifi)^ f (expcc>^)Td(Af^J. (27) 

The right hand side is the Hirzebruch-Riemann-Roch number of M^. The only factor in the integral which 
depends on /i is the symplectic form, everything else depends only on the region containing fi. Thus by 
l)25|l m\{p) is a polynomial function of ^ on A^ as with the Duistermaat-Heckman measure. I 
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Remark 3.3 This proof implies that the optimal domains of polynomiality for the multiplicity function 
must be unions of the Ai's. Guillemin, Lerman and Sternberg ^14-! have shown that this partition is opti- 
mal for the Duistermaat-Heckman measure by showing that the difference between the polynomials in two 
adjacent regions is nonzero. We conjecture that this partition is also optimal for the multiplicity function. 
This has been confirmed up to SL4C. 



As further evidence for the conjecture, we note that on a given domain, the weight polynomial and the 
Duistermaat-Heckman polynomial in H2t)|) have the same leading term since 

d 

Td(Af^) = l + ^r, (28) 
with Tj e 57^-' (Af^(j) in the de Rham complex. 

Remark 3.4 It is a very interesting open problem to count the regions in the permutahedron subdivided 
according to Theorem This is the analog of Kirillov's question for the Kostant partition function 

mentioned in Section \l.4\ We have determined all the region counts for SL4C in Figure\^ 



There are many links between the weight multiplicities and the Duistermaat-Heckman function. For ex- 
ample, Dooley, Repka and Wildberger [5] provide a way to go from the weight diagram for A to the 
Duistermaat-Heckman measure for 0\+s' 

f^?s= E -a(/3)/F"- (29) 

13 weight of Vx 

Also, if f is the Lebesgue measure on t*, normalized so that the parallelepiped given by the simple root 
vectors has unit measure, we define the Duistermaat-Heckman measure to be the product f^^i'. Now for 
each n G N construct the discrete measure 

^" ^ dimV„A ^ rnnx{P)Sp/n (30) 

13 weight of Vnx 

where Sx is a point mass at x and VnX is the irreducible representation of su{k) with highest weight nX. 
Then Heckman |16| proved that Vn converges weakly to the Duistermaat-Heckman measure as n — > cxd. 

Furthermore, the Duistermaat-Heckman function from above can be computed in the following way (see 
|14p: its value at a point (A, /3) is obtained by using Kostant 's multiplicity formula with 6 = and a 



deformation of Kostant's partition function that takes the volume of the polytopes {{ka)aeA+ G 
SasA ^aOi = v} instead of their number of integral points. 



|A+| 
^>0 



4 The Kostant arrangements 

In this section, we will construct a hyperplane arrangement whose regions are also domains of polynomiality 
for the multiplicity function. This partition into domains will be unlike the ones obtained in Remark 12.51 
and Theorem 13.11 in that it is not invariant under rescaling A and /?. We will deduce the form of this 
arrangement from a closer look at Kostant's multiplicity formula © and its chamber complex defined in 
Section O 

Lemma 4.1 The set of normals to the facets of the maximal cones of the chamber complex of the Kostant 
partition function of An consists of all the conjugates of the fundamental weights. 

Proof. The facets of the maximal cones of the chamber complex span the same hyperplanes as the facets 
of the base cones whose common refinement is the chamber complex. Base cones correspond to sets of n 
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linearly independent positive roots. Fixing a particular base cone spanned by {71, ...,"/„}, consider the 
undirected graph G on {1, . . . , n + 1} where is an edge if — Cj = 7™ for some m. The fact that the 
7j's are linearly independent implies that G has no cycles. So G is a forest, and since it has n + 1 vertices 
and n edges (one for each 7^), it is actually a tree. Suppose now we remove 7^ — Cg — et and want to find 
the normal of the hyperplane spanned by the other 7,j's. The graph G with the edge (s, t) removed consists 
of two trees Ti and T2. List {l,...,n+l} in the form 

a : ii, 12, . . . , ij-i, s, t, . . . , in+1-2 
■> , ' , ■■ 

vertices of Ti vertices of T2 

where we will think of ct as a permutation in one-line form. 

Now let a'j^ = ea-(i) — eo-(i+i) and note that a'j ^ Cs — et = jj. The set {a[, . . . , a'^} is a root system basis 
because it is the image under the action of of the original simple roots ai — Ci ~ e^+i. Observe that 
every edge in Ti can be expressed as a sum of a'j^, . . . , a'j-i, and every edge in T2 as a sum of ttj+i, • ■ • , a'm 

so that all 7^ 's in {71 , . . . , 7^ , . . . , 7„ } can be expressed as linear combinations of , . . . , , . . . , . The 
normal for the corresponding hyperplane will therefore be the jth fundamental weight lu'j for the basis 
{a[, . . . , a'„} = a • {ai . . . , an}- 

Conversely, given any fundamental weight loj for the root system basis cr • {ai . . . ,q;„} (or equivalently, 
(T~^ • ujj, where ujj is the jth fundamental weight for the standard simple roots), we want to show it can 
occur as the normal to a hyperplane. Let iJ be a hyperplane separating the standard positive roots from 
the negative ones. For each a[ — a ■ ai, we can pick a sign such that e^a^ is on the positive side of H. 
Hence {eia'i, . . . , Sna'^^} is a linearly independent subset of the set of standard positive roots, and thus it 
corresponds to one of the base cones of . The corresponding graph is a path since we have a system of 
simple roots (up to sign reversal). Removing Eja'j and applying the above procedure with the order given 
by the path gives that loj occurs as the normal of the corresponding hyperplane. I 

To compute multiplicities for 5[fcC using Kostant's formula, we look at the points (t(A + S) — {P + S), as a 
ranges over the Weyl group 6k- Some of these points will lie inside the chamber complex for the Kostant 
partition function and we compute the multiplicity by finding which cells contain them and evaluating the 
corresponding polynomials at those points. Starting with generic A and (3, none of the points (j{X+d) — {(3+6) 
will lie on a wall of the chamber complex of the Kostant partition function, and if we move A and /? around 
a little in such a way that none of the a{X + 6) — {(3 + 6) crosses a wall, we will obtain the multiplicity for 
the new A and (3 by evaluating the same polynomials. So there is a neighborhood of (A, /3) on which the 
multiplicity function is given by the same polynomial in variables A and (3- 

Lemma 14.11 describes the walls of the chamber complex for the Kostant partition function in terms of the 
normals to the hyperplanes (though the origin) supporting the facets of the maximal cells. Now a point 
(t(A + 5) — {f3 + S) will be on one of those walls (hyperplane though the origin) when its scalar product with 
the hyperplane's normal, say 0{u!j), vanishes, that is when 

(a{\ + S)-{f3 + 5),e{LJj)) =0 (31) 

For any A, consider the arrangement of all such hyperplanes for I < j < k and a, 9 &k- For (3 and /?' in 
the same region of this arrangement and any fixed a G &k, the points cr(A + (5) — {(3+6) and a{X+6) — {(3' + 6) 
lie on the same side of every wall of the chamber complex for the Kostant partition function. Figure 13 (on 
the left) shows the arrangement we get for A — (11, —3, —8) in A2, with and without the weight diagram. 

In view of the invariance of the multiplicities under the action of the Weyl group, Kostant's formula has 
to give the same thing if we replace (3 by '0(/3)i ^ G &k- Replacing (3 by ip{l3) in the equations 1)31(1 above 
yields another hyperplane arrangement, which we will denote by A^^'^ - Hence for each A, we get a family 
of arrangements indexed by t/i G &k, which we will call the Kostant arrangements for A. Figure 01 (on the 
right) shows the superposition of those arrangements for A = (11, —3, —8) as ^p ranges over the Weyl group. 

Suppose the chamber complex for the Kostant partition function for slfeC has r{k) full dimensional cones. 
We will choose a labeling of these regions with the integers 1, . . . , r{k) once and for all, and let the associated 
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polynomials be pi , . . . ,Pr{k)- Recall, these are polynomials of degree 2 ^) "-"^ ^^"^ subspace xi + - ■ ■ + Xk — 
of MJ'. We will also label the exterior of the chamber complex by and let its polynomial be the zero 
polynomial pq . 

Definition 4.2 For generic A and (3, let v^\\,j3) (or just ) be the label of the region containing the 
point a{X + S) ~ (V'(/5) + ^) (this label is unique for generic A and j3). Define the type of A and f3 to be the 
vector 

Type('^)(A,/3) = («W).ee., 
for some fixed total order on 6fc. Furthermore, define 

Px\P) = E (-l)''^V.<.) + - (32) 

Lemma 4.3 p'^^ is a polynomial function on the interior of the regions of A^^^ and coincides with the 
multiplicity function there. 

Proof. For fixed A, the type of points along a path between two /3's in the interior of the same region 
of will remain the same by definition of the Kostant arrangement (because no (t(A + 5) — (i^iP) + S) 
crosses a wall along that path). ■ 

The reason why Lemma 14.31 is restricted to the interior of the regions is that while polynomials for adjacent 
regions of the chamber complex for the Kostant partition function have to coincide on the intersection of 
their closures, there is a discontinuous jump in the value of the Kostant partition function (as a piecewise 
polynomial function) when going from a region on the boundary of the complex to region (outside the 
complex) . 

Remark 4.4 Given a rational polytope Q of dimension d in R" and t E N, denote by tQ the polytope 
obtained by scaling Q by a factor of t. Ehrhart showed that the function t 1-^ \tQ n Z"|, counting the 
number of integer points in tQ as a function oft, is a quasipolynomial of degree d, and a polynomial of degree 
d if Q is integral. This function is called the Ehrhart (quasi) polynomial of the polytope Q. Furthermore, 
the leading coefficient of the Ehrhart quasipolynomial is the d-dimensional volume of Q. It can be shown 
that for every fixed A and (3, and any ip, the function t 1— > Type'''^'' {t\, tj3) in the nonnegative integer variable 



4 The Kostant arrangements 



15 



t eventually stabilizes as t grows (in a way that depends only on k and not on A and j3). This can be used 
to give a proof that the Ehrhart functions of the Gelfand-Tsetlin polytopes GTa,/3 are polynomial (we omit 
the proof since we prove something stronger in Corollary \5.S\ below). 

In the definition of the Kostant arrangements above, a lot of the hyperplanes are redundant. We simpHfy 
here the description of these arrangements. Since everything occurs on the subspace Xi + - ■ ■ + Xk = of R*"', 
so that (cr(A + (5) — + (1, 1, ■ ■ • , 1)) = 0, we can regard the normals to the hyperplanes up to adding 
multiples of (1,1,.. .,1) without changing the Kostant arrangements. So we can use ujj = cl'j + |:(1, 1, . . . , 1): 

• + ej , 




which is more convenient than ujj for what follows. The hyperplanes of Ay^ then have the form 



- (a(A + ,5)-(^(/3) + <5),0(cD,)) 

= ((T(A)-VX/3)+(T(<5)-(5,ee(i) + --- + ee(j)) 

= ((Acr-i(i) - + <5cr-i(j) - (5j)i=i,...,fc: H ^ ^e{j)) 

j 



i=l 



H ^ Pip-HO{j)) — + h ACT-i(e(j)) + ^((5a-i(e(i)) - <^e(i)) • (33) 



At this point we can get rid of the permutations since only the subsets 0{{1, . . . ,j}), ijj^^9{{l, . . . , j}) 
and a~^9{{l, . . . , j}) are important and not the order of their elements. They can be any subsets since 
ip~^9, a~^9 and 9 can be any three permutations of &k- Because the (3i, the Ai and the 6i sum up to 
zero, replacing these subsets by their complements gives the same hyperplane. This proves the following 
proposition. 

Proposition 4.5 The hyperplanes of the Kostant arrangements are defined by the equations 

j 

Pu, + --- + Pu,=K+--- + + -5^,), (34) 

i=l 

where U = {ui,...,Uj}, V = {vi, . . . ,Vj} and W = {wi,...,Wj} range over all j- element subsets of 
{l,...,k} and j < lk/2\. 

We will call the correction term involving only S in a hyperplane given as in H34|) , the 5 -shift: 

3 

shm{V,W)=Y,{Sv,-S^^). (35) 

1=1 

Remark 4.6 For fixed U , we get a series of parallel hyperplanes, and we can determine which are the 
outer ones because they correspond to maximal and minimal sums of Ag. . Since the Ai > A2 > • • • > Afe, 
they are 

(im+--- + Pu,= Ai + --- + A, + shift({l,...,j},M^) 

(36) 

f3u,+--- + Pu,= Afe_j+i + ■ • • + Afe + shift({fc- j + l,...,fc},W^). 

Note that since the coordinates of 5 are decreasing, shift({l, . . . , j}, W) > and shift({fc— j+1, . . . , fc}, W) < 
for all W. 
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5 POLYNOMIALITY IN THE CHAMBER COMPLEX 



We conclude this section by relating the domains given by the Kostant arrangements and those given 
by Theorem 13.21 by showing that the hyperplanes supporting the facets of the domains are precisely the 
hyperplanes of the Kostant arrangements without the (5-shift factors. 

Proposition 4.7 The supporting hyperplanes of the facets of the top- dimensional domains of the permu- 
tahedron for generic X are the hyperplanes 

Pu,+--- + f3u,=K + --- + , (37) 

for 1 < j < [fc/2j and U = {ui, . . . , uj}, V = {vi, . . . , Vj} ranging over all pairs of j-element subsets of 
{1,...,A:}. 

Proof. Theorem 13.11 gives the walls supporting the facets as the convex hulls of W ■ <j{X), where cr(A) is a 
point of the Weyl orbit of A, and is a parabolic subgroup of the Weyl group. For 6k, those subgroups 
permute two complementary sets of indices independently. If U is one of those sets of indices, with \U\ = j, 
and A^j , . . . , A„j the coordinates of o'(A) in those positions, then the hyperplane supporting W ■ cr(A) is 

+ • ■ • + = A„, + • • • + A„^- . (38) 

Had we chosen the complement of U instead with the remaining A^'s, we would have gotten the same 
hyperplane in the subspace xi + • ■ • + xj, = of M*^ since the A^'s and the Pi's sum up to zero. ■ 

We can obtain the following corollary without using the full description of the domains of the permutahedron 
obtained by symplectic geometry means in Theorem 13. II 

Corollary 4.8 The hyperplanes supporting the facets of the permutahedron for a generic A are 

Pui + ^ (3uj = Ai H \- Xj , . 

+ • ■ • + Pu, = Xk-j+i + --- + Xk ^ ' 

for 1 < j < [k/2\ and U — {wi, . . . , Uj} ranging over all j-element subsets of {1, . . . , k}. 

Proof. We remark that the "shell" of the weight diagram is just a permutahedron, whose facets can easily 
be described in terms of permutations (see O p. 18]). For U C {1, . . . , fc}, construct a fc-vector by putting 
the first \U\ Ai's in the positions indexed by U and filling the other positions with the remaining elements. 
Then act by the subgroup of 6^ that permutes the elements in positions U and {1, . . . , k}\U independently 
to get a facet as the convex hull of the points of this orbit. The affine span of this facet is the hyperplane 

+ ■ ■ • + - Ai + • • ■ + A^- . 

By choosing the last \U\ Ai's instead, we get the hyperplane supporting the opposite parallel facet. These 
are the outer hyperplanes from (|36|l without the shifts. Remark 14.61 also implies these outer hyperplanes 
actually lie outside the permutahedron. ■ 

Remark 4.9 The Weyl orbits ofujj andujk-j A < j < A; — IJ/orsI^C (type Ak-i) determine the same set 
of directions, since LUk-j is —ujj with the coordinates in reverse order. So the Weyl orbits of loi, . . . ,L0[k/2\ 
already determine all the possible normals to facets of the permutahedron ( and the hyperplanes of the Kostant 
arrangement) . 



5 Polynomiality in the chamber complex 



Theorem 12.11 allowed us to write the multiplicity function as a partition function, which is therefore 
quasipolynomial over the convex polyhedral cones of the chamber complex C^'^-'. On the other hand. 
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for each dominant weight A, Theorem 13.21 shows that the partition of the permutahedron from Theorem l3.1l 
gives domains over which the muhiphcity function is polynomial in [3. We show here that the quasipolyno- 
mials attached to the complex C^*"' are actually polynomials, so that the multiplicity function is polynomial 
in both A and [3 over the cones of the complex. 

The union of the cones of the complex C^^'^ is the cone 

r^'^ - U ^ conv(efc • A) , (40) 

AG Co 

where Co is the fundamental Weyl chamber. 

We can lift the partition of the permutahedron from Theorem 13. II to (A,/3)-space by lifting the wall 

conv(M^ ■ (t(A)) 

to 

y {A} X conv(W^ • cr(A)) . 

AG Co 

This gives a partition T'^^'^ of the cone T'^'^' into convex polyhedral cones, and Theorem 12 . II implies that the 
multiplicity function is quasipolynomial over the cones of rW. We recover the domains from Theorem 13. II 
by intersecting T^'') with L(A) from equation H19|l . Our reason for introducting T*^'') rather that working 
with the complex C^*^^ is that Proposition 14. 71 lets us describe the hyperplanes supporting the facets of the 
cones of T'^*^-' easily. Indeed, if 

/3ui + • • ■ + Pu, ^\v^+■■■ + Xv, 

supports a wall conv(T4^ • o'(A)) for fixed A, then the wall UagCo'I'^J' ^ conv(T4^ • o'(A)) is supported by the 
hyperplane 

/3«i + • • • + = A^j + • • • + 
in (A, /3)-space, where we now think of A as variable, just like /?. 

The last tool we need is a lifted version of the Kostant arrangements. Recall from (|33|l that the Kostant 
arrangement A'"^^ has the hyperplanes 

j 

Pi:-^(e(i)) + ^ Pi:-^(0(3)) = ^a-^(e(i)) + \- ^c-^(e(])) + '^^{^<y-^(e(i)) - ^9(1)) , 

1=1 

as 6 ranges over 6^. We will denote by A^^"^ the arrangement with hyperplanes 

j 

Pi^-HsW) + ^ (^i^-HHj)) = ^<y-Heii)) + ^ ^<y-H9U)) + X^^'^'^^M^W) " ^e{i)) , 

1=1 

where we now think of A as variable, and 9 ranges over &k as before. The definition of Pj''''' (equation H32I) ') 
and Lemma |4.3I generalize to give us a piecewise polynomial function P^"^) in A and P that expresses the 
multiplicity function as a polynomial on the interior of the regions of the arrangement A^"^^ . 

Theorem 5.1 The quasipolynomials determining the multiplicity function in the cones of T^^^ and C^*^-* 
are polynomials of degree C^j^) in the j3i, with coefficients of degree ('^2^) in the Xj. 

Proof. We will show that for each cone C of T'--''^ we can find a region R of the Kostant arrangement 
^(0) (for any -0), such that C d R contains an arbitrarily large ball. Then p('^) and the quasipolynomial 
in C agree on the points (A, /3) in that ball for which (A, (3) G Aw x Aw and A — /? G Ar (the points 
corresponding to allowable pairs of a dominant weight and a weight of its irreducible representation) . The 
quasipolynomial must therefore be polynomial on those points. The degree bounds follow from r!,emark ll.7l 
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6 Factorizations of weight polynomials 



By the remarks preceding this theorem, the hyperplanes supporting the facets of the cones of T'*^^ are 
exactly the same as the hyperplanes of the Kostant arrangement A'^'^'^ with the shifts removed. If we 
deform A^^^ continuously to make the shifts zero (by multipliying them by t and letting t going from 1 to 
0, for example), the final deformed arrangement is a partition of T'^'^^ that refines T^''\ Let R be any region 
of A'^^'^ whose deformed final version is contained in C. Consider a ball of radius r inside the deformed 
image of R, and suppose it is centered at the point x. If s is the maximal amount by which the hyperplanes 
of the Kostant arrangement are shifted, then R contains the ball of radius r — s centered at x, and so does 
C n R. Since C is a cone, we can make r arbitrary large and the result follows since s is bounded. 

We get the same result for the complex C''^^ by passing to its common refinement with T''^^. ■ 

Recall from Section |21 that the weight multiplicity function and the Duistermaat-Heckman function have 
the same leading term. In particular, the degree of the multiplicity function is at most the upper bound 
on the degree of the Duistermaat-Heckman function. For a torus T acting on a symplectic manifold M, 
the latter is known to be (dimM)/2 — dimT. In our case, M is the coadjoint orbit 0\ and dimT = fc — 1 
since T is the set oi k x k traceless diagonal Hermitian matrices. The dimension of 0\ is k^ ~ k ^ k{k — 1) 
for generic A, but for nongeneric A, we can get more precise bounds on the degrees. Since the coordinates 
of A are decreasing, it has the form 

/ci times k2 times 

where > 1^2 > ■ ■ ■ > i^i and the kj sum up to k. In this case, one can show that dimO;^ = ^ ' 
that the weight multiplicity function for that A is piecewise polynomial of degree at most 

+ (41) 

For s^C, for example, we get at most cubic polynomials for generic A, at most quadratic polynomials for 
A with exactly two equal coordinates, at most linear polynomials for A with two pairs of equal coordinates 
(A of the form (:/, v, —v, —v)) and constant polynomials for A with three equal coordinates (A of the form 
(3i/, — i/, —h', —v) or [v, V, — 3i^)). 

We can also deduce from Theorems 12 . II and 15.11 that the multiplicity function for type A exhibits a scaling 
property in the following sense. 

Corollary 5.2 Let T he the set {(A, /?) £ : A — /3 G A/?}. For any generic (A, /?) £ T, we can find a 
neighborhood U of that point over which the function 

(A,/3,i) e (C/nT) xN^mtA(t/3) (42) 
is polynomial of degree at most 2('^2^) in t and ('^j^) in the A and (3 coordinates. 

Proof. Let (A,/3) e T. For U sufficiently small, the points {(tA,t/3) : i G N} lie in the same cone of the 
chamber complex C'-'^-', and for t S N, <A and tj3 are points on the weight lattice with their difference on 
the root lattice. Hence the corresponding multiplicities are obtained by evaluating the same polynomial at 
those points. ■ 

Remark 5.3 This corollary implies in particular that the Ehrhart functions (see Remark \4.4^ of the 
Gelfand-Tsetlin polytopes GT^.^ are always polynomial, even though the polytopes are not always integral 
(see f^;. 



6 Factorizations of weight polynomials 

In this section we use the explicit relation between the hyperplanes of the Kostant arrangements and the 
supporting hyperplanes of the partitioned permutahedron to identify certain factors in the weight polyno- 
mials. As mentioned in the introduction, Szenes and Vergne |31j have recently observed this factorization 
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phenomenon for general partition functions. The quasipolynomials associated to the partition function's 
chamber complex exhibit a certain number of linear factors that vanish on hyperplanes parallel and close 
to those supporting the walls of the complex. In our case it is unclear how to deduce the form of the 
walls of the complex C^''^ from the complex of the partition function given by matrix E/^ in Section |21 
We are however able to deduce similar results from the Kostant arrangements and the description of the 
hyperplanes supporting the walls partitioning the permutahedron from Section 0] 

6.1 On the boundary of the permutahedron 

We have seen in Proposition that each facet of the permutahedron is parallel and close to a hyperplane 
of a Kostant arrangement. This means that the domains of polynomiality of the weight diagram that are 
on the boundary of the permutahedron overlap with regions of the Kostant arrangement, but can't coincide 
because of the shifts caused by S. We can use this to our advantage to show that those weight polynomials 
have to factor somewhat. The reason is that two polynomials give the weights in the overlap: the one 
attached to a cone of the chamber complex obtained from writing the multiplicity function as a single 
partition function, and one, P^'^\ coming from Kostant's multiplicity formula. Because the overlap isn't 
perfect, the polynomial from Kostant's formula is valid on a region that goes outside the weight diagram 
and must therefore vanish there. The purpose of this section is to make precise this phenomenon and 
quantify it. 

Definition 6.1 For fixed X, consider the hyperplane 

H : (3ui + • • • + l^uj = Ki + ■ ■ ■ + Kj 

where 1 < j < [fc/2j and U — {ui, . . . , Uj}, V — {vi, . . . , Vj} are j-element subsets of {1, . . . , k}. We will 
call the polynomial 

7c/,y(A) ^I3u,+---+ (3u, ~Xv, A„^. e Z[f3] (43) 

the defining equation of H . For variable A, we also define 

lu,v = + ■ ■ • + Pu, - K e P] ■ (44) 

Theorem 6.2 Let R be a domain of polynomiality for the weight diagram of the irreducible representation 
of slfeC with highest weight X, and pn he its weight polynomial. Suppose that R has a facet lying on the 
boundary of the permutahedron for X that has Oiujj) as its normal vector, for some 9 € 6^. If "f = 7(A) is 
the defining equation of the hyperplane supporting that facet, then p^ is divisible by the j{k ~ j) — 1 linear 
factors 7 + 1, 7 + 2, . . 7 + j(/c - j) - 1, or 7 - 1, 7 - 2, . . 'y~j{k-j) + 1. 

Observe that this is invariant under replacing j by k—j, which is a consequence of the remark in Rcmark l4.9l 
By that remark, we can therefore restrict ourselves to 1 < j < [fc/2j. 

Proof. Suppose the hyperplane supporting the facet F of R on the boundary of the permutahedron has 
normal 0{ujj). By Remark l4.(j| and Proposition this hyperplane is either 

Peii) + ■ ■ ■ + [36(3) = Ai + --- + Aj or (45) 
Peii) + ■ ■ ■ + f^eu) = Afe_j+i + --- + Afe. (46) 

Suppose it's the first one (the argument is the same for the second one). From Proposition 14.51 we know 
that in the Kostant arrangements, we have the hyperplanes 

Pe(i) + ■■■+ I3eu) - Ai + ■ • • + A, + shift({l, . . . , j}, W) , (47) 

for W ranging over j-element subsets of {1, . . . , /c}. We want to identify a region R' of one of the Kostant 
arrangements that overlaps with R and extends beyond the boundary of the weight diagram as far as 
possible. Note that although the exterior walls of R and R' have to be parallel, the interior walls don't. 
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6.1 On the boundary of the permutahedron 




First we find a hyperplane of one of the Kostant arrangements parallel to F and outside the permutahedron. 
Recall from Remark 14.61 that hyperplanes of the form (|47|l have nonnegative 5-shifts. For positive (5-shift, 
the corresponding hyperplane lies outside the permutahedron. In fact, we would like to maximize 

3 3 

shift({l, ...J},W)^Y.^^-Y. (48) 

i=l 4=1 

because this will determine how much pR factorizes. The first sum is as large as possible because it is the 
sum of the first j coordinates of 5 — |(fc — 1, fc — 3, . . . , — (fc — 3), — (fc — 1)). Since j < [fc/2j , we can pick 
W disjoint from {1, . . . ,j}. Picking W = {k — j + 1, . . . , k} means the second sum consists of the last 
(and smallest) entries of 6. Thus (fc — l)/2, (fc — 3)/2, . . ., (fc — 2j + l)/2 appear in the first sum and their 
opposites in the second. The maximal shift is then 

shiftC- )(,) = 2 + ^ + • • • + ^^f^) = - j) . (49) 

Suppose that H(X) is the hyperplane with this maximal shift (at distance j{k^j) outside the permutahedron 
and parallel to F) and that it belongs to the Kostant arrangement A^^\ 

The second step is to find a region R' of -4^''''' with a facet on H{X) that overlaps with R. If we replace A 
by a multiple mX of itself, the partition of the permutahedron simply scales up by a factor of m, and the 
polynomials attached to the regions, as polynomials in A and (3, remain the same (because the cells of the 
chamber complex C''^-' are cones). The hyperplanes of the Kostant arrangements almost scale, except for 
the (5-shift factor. Those shifts preserve the distance between the hyperplanes and the ones supporting the 
facets of the permutahedron, even as the regions grow since the separation between parallel hyperplanes of 
increases. Hence for a large enough multiple of A, one of the regions R' of Al^f] with a facet on iJ(mA) 
will overlap with mR. From now on we'll assume that A has been replaced by a suitably large multiple of 
itself. 

We are now in the setup of the above picture. The polynomials in R and P^'^) on R' both give the 
multiplicities in the interior of their respective regions, and hence they are equal provided that R Cl R' 
contains sufficiently many points. We can assume that we have scaled A sufficiently above so that this is 
the case. Since P'^^'^ has to vanish outside of the permutahedron, it will vanish on the intersection of R' 
with the hyperplanes 

Pe{i) + ■■■+ /3eo) = Ai + ■ • • + Aj + 1 

f3g(i) + ---+ Peu) = Ai + . • • + Aj + shift^'"^'^) (j) - 1 . 

If the intersection of R' with these hyperplanes contains sufficiently many points (again, we can scale A so 
that this is the case), P^''^' will have the defining equations of those hyperplanes as factors, and hence so 
wiU PR. 
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Here we have assumed F is defined by and 7 = ^u,v{^) for U ~ 6'({1, . . . , j}), V — {!,..., j}. 
If F is defined by (|46|) . we get the same C/ but = {fc — + 1, . . . ,/c}, and the defining equations 
7,7-l,---,7-shift('"^'^)(j) + l. ■ 

We can lift this result to the weight polynomials associated to the cones of the chamber complex C^*'-'. This 
will allow us to think of the linear factors dividing the weight polynomials as polynomials both in A and (3. 

Corollary 6.3 Let r be the cone of C'^^^ whose intersection with L{X) gives domain R in the previous 
theorem, and pr its associated weight polynomial. If ^u.vW + c divides pR, then ^u,v + c divides pr- 

We will call these families of linear factors, parallel linear factors. This shows that the smallest number of 
parallel linear factors is obtained when considering facets of the permutahedron normal to a permutation 
of LOi or LOk-i- In this case, we get k — 2 factors in the weight polynomials of the boundary regions on those 
facets. For j = [fc/2j, we get a maximum of [■|J(fc — [|J) — 1 ^ parallel linear factors. Since the pr 
have degree at most C^j^) ^ fc^/2 in /? (regarding A as a parameter), we get linear factors accounting for 
about half the degree of the weight polynomials for those facets. 

The fundamental weight LUj = ei + 62 + ■ ■ ■ + ej ~ |(1, !,...,!) has an orbit of size (^), and thus there 
are that many facets having a permutation of luj as an outer normal (the opposite parallel facets have the 
permutations ofcuk-j as normals). So the permutahedron for generic A = (Ai, . . . , A^) has X]j=i (j) — 2''' — 2 
facets, most of which have normals corresponding to central values of j (i.e. close to [fc/2j). The following 
table gives the minimum numbers of parallel linear factors for different values of k and j. In parentheses are 
the numbers of facets having a permutation of ujj or tOk-j as a normal. For example, in slsC, the maximal 
degree of the weight polynomials is 21 and we expect that the polynomials of regions with a facet on any 
of 112 of the 254 facets of the permutahedra to have 14 parallel linear factors. 





#(facets) 


maxdeg(pH) 


J-1 


J-2 


J - 3 


J = 4 


skC (A2) 


6 


1 


1 (6) 








skC (A3) 


14 


3 


2 (8) 


3 (6) 






SI5C (Ai) 


30 


6 


3 (10) 


5 (20) 






sleC (As) 


62 


10 


4(12) 


7 (30) 


8 (20) 




sItC (Ae) 


126 


15 


5 (14) 


9 (42) 


11 (70) 




skC (Ar) 


254 


21 


6 (16) 


11 (56) 


14 (112) 


15 (70) 


SigC (As) 


510 


28 


7(18) 


13 (72) 


17 (168) 


19 (252) 



Theorem 16 . 21 onlv depends on the fact that using our description in the previous section of the walls of the 
chamber complex for the Kostant partition function and the combinatorial description of the permutahedron 
(Lemma I4.1|l , we can argue that there will always be hyperplanes of Kostant arrangements parallel and 
close to the facets of the permutahedron. In order to extend the factorization phenomenon inside the 
permutahedron, we will need to use the complete description of the domains of polynomiality for the 
weight multiplicity function, obtained by symplectic geometry means in Theorem 13. II 

6.2 Inside the permutahedron 

We already discussed at the end of the previous section that the hyperplanes supporting the walls parti- 
tioning the permutahedron are precisely the hyperplanes of the Kostant arrangements without the shift 
factors. We will take advantage here of overlaps between the improved domains of Section Inland regions of 
the Kostant arrangements to show, not that the weight polynomials themselves factor, but rather that as 
we jump between two adjacent domains, the difference in the corresponding weight polynomials exhibits 
parallel linear factors. Given a facet between two adjacent domains of the permutahedron, we will see that 
we are able to find two hyperplanes of Kostant arrangements parallel to it and at maximal distance on 
either side of it, and deduce from this a number of parallel linear factors of the polynomial jump. 
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6.2 Inside the permutahedron 



Definition 6.4 We will say that two domains are adjacent if they have the same dimension and a facet of 
one is a subset of a facet of the other, or equivalently if they intersect in a nonempty polytope of dimension 
one less. 



Theorem 6.5 Let Pi and P2 be two adjacent full dimensional domains of polynomiality of the permuta- 
hedron for a generic dominant weight A of slkC, and suppose that the normal to their touching facets is in 
the direction cr(iOj) for some a G 6fe. If Pi and P2 are the weight polynomials of Pi and P2, and 7 is the 
defining equation of the wall separating them, then the jump pi — p2 either vanishes or has the j{k — j) — 1 
linear factors 

(7 - s- + 1), (7 - s- + 2), . . . , 7, . . . , (7 + s+ - 2), (7 + s+ - 1) 
for some integers s~ , s^ > satisfying 

s- +S+ =j{k~j). (50) 



Proof. Suppose the touching facets of Pi and P2 he on the hyperplane 

/3«i + • • • + = + • • • + . 

Then among the Kostant arrangement hyperplanes 

+ • ■ • + = A„, + • ■ • + A„^ + shift(F, W), WC{l^...,k}, \W\ = j , 

we can find a pair for which the ^-shift is minimal and maximal by picking appropriate subsets W. Clearly, 
the minimal shift — will be nonpositive, and the maximal shift s+, nonnegative. In fact, 



s+ 



s 



so that 

s+ + s 



since 6 — — (5icvcrsc_ p^^. ^ odd, the Si are integral, and hence so are s" and s+. When k is even, the Si 
are half-integers with odd numerators. Since we are adding/subtracting an even number of them {2j) to 
compute the shifts, we again get that s^ and s+ are integers. 

We can find regions Qi and Q2, Ri and R2 of Kostant arrangements as in the following diagram. We 
will think of these regions as open convex polytopes because in Lemma 14.31 the polynomials giving the 
multiplicities on the regions of the Kostant arrangements are only valid in the interior of the regions. 



i=l i=l J i=l i=l 

(j j \ j 3 

1=1 1=1 / i=l 1=1 



max > 5w — min > Sw = 2max > Sw = i(k — j) (51) 



w w w 

i=l 2 = 1 2 = 1 



7 The chamber complexes for SI3C {A2) and s^C (A3) 
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We will let the corresponding polynomials, as given by P^"^") in equation H32|) . be qi and 92, ''i and r2 
respectively. Since wc can assume that we have scaled A sufficiently (as in Theorem 16. 2|) . we have that 
qi = Pi = fi and 52 — P2 = ^2, since Qi H Pi n i?i and (52 H P2 H R2 are large. Furthermore, pi and 52 
agree on Pi n Q2 and similarly, p2 and ri agree on P2 fl i?i . Since Pi n (32 and P2 n Pi contain enough 
lattice points on the bounded hyperplanes (dotted lines in the diagram) , the differences pi — (72 and p2 — ri 
have to vanish on those hyperplanes. Hence 

Pi -92 = (7-s- + l)(7-s- + 2)---(7-l)-/ii 

(52) 

P2-ri = (7+l)(7 + 2)---(7 + s+-l)-/i2 

for some polynomials hi and ft-2, unless pi = 92 or p2 = hi which case pi — P2, since pi — ri and p2 = 92- 
If we assume that pi P2, we have that 

P1-P2 = (7-s- + l)(7r + 2)---(7-l)-'ii , , 

(53) 

P2-P1 - (7+l)(7 + 2)---(7 + s+-l)-/i2 

and since pi and p2 have to agree on the lattice points on the wall between Pi and P2, their difference is 
also divisible by 7. Hence we get 

Pi - P2 = (7 - s")(7 - + 1) • • • 7 • • • (7 + s+ - 1)(7 + s+) • /13 (54) 

for some h^. M 

Remark 6.6 As in Corollary 1 6'. 51 we can lift this result to the weight polynomials Pr associated to the 
cones of the chamber complex C*^*^-* and regard the parallel linear factors as polynomials in both A and (3. 



7 The chamber complexes for ^[sC (^2) and s^C (A3) 

In this section we explicitly compute the chamber complexes for fc = 3 and k — 4. For A; = 4, we find 
that the chamber complex does not optimally partition the domains of polynomiality for the multiplicity 
function. In Theorem l7.4l we prove that the optimal glued complex does agree with Theorem 13 . 21 for k — A. 



7.1 The chamber complex for sl^C (A2) 



Using the procedure described in Section |21 to write down the multiplicity function as a single partition 
function in the case k — 3 {A2) gives that 

mA(/3)=0B3 (ssQ)) (55) 
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7.1 The chamber complex for bIsC {A2) 



with 



/I 1 \ 
0-11000 
10 10 
0-10010 

\ -1 1 / 



and Bs 



/ Ai - A2 \ 

2X2 -Pi 

/3i + /32 + Ai 
A2-/3i 

V A2 - /32 / 



(56) 



We can compute the chamber complex associated to E3 and intersect it with the space 



Sal J : AGM^/3eM^Al+A2 + A3 = 0,/3l+/32+/33=0i> . 



(57) 



In that space we can apply 



b; 





1 


6 


2 


3 


1 


1 \ 


1 




-3 


2 


3 


1 


1 


9 




-3 


-1 


3 


-5 


4 






-3 


-1 


3 


4 


-5 / 



to rectify the cones of that complex to obtain C^^\ The full dimensional cones of C'^^^ are given by 







T3 


= pos(6, oi, 


C2 


C3) 


T6 


= pos(6, ai 


0-2,03) 


n = pos(6, 01,02,03) 


T2 = P0S(6, Cl,C2,C3) 


Ti 


= pos(6, 02 


Cl 


C3) 


T7 


— pos(6, ai 


as, C2) 






T5 


= pos(6, 03 


Cl 


C2) 


T8 


= pos(6, 02 


a3,ci) 



where the rays are 

ai = [2, 
a2 = [2, 
as = [2, 



-1,-1,2,-1,-1] 
-1,-1,-1,2,-1] 
-1,-1,-1,-1,2] 



6= [1,0,-1,0,0,0] 



ci = [l,l, 

C2 = [1, 1, 
C3 = [1,1, 



-2,-2,1,1] 
-2,1,-2,1] 
-2,1,1,-2] 



(58) 



(59) 



(60) 



The cones above are grouped into orbits under the action of the symmetric group S3 on the /3-coordinates. 
In general the set of cones of won't be closed under the action of Sfe on the /3-coordinates, even though 
the multiplicities under the &k should be invariant. 

We can get the polynomial Pi corresponding to r, easily through interpolation, using for example the 
Kostant partition function for A2 which has the simple form 



K{a, —a + b. 



min{o, 6} 




1 



if o, 6 G N , 
otherwise. 







P3 = 


l + Ai 




P6 


= 1 + Al -f A2 + As 


Pl = 1 + A2 - A3 


P2 = 1 -1- Al - A2 


Pi = 


1 + Ai 




P7 


= 1 + Al + A2 -f /32 






P5 = 


1 + Ai 


-/?3 


P8 


= 1 + Al + A2 -h /3i 



Note that even though they highlight the symmetries in the /3i's, those polynomials are a little ambiguous 
since they are defined up to the relations Ai + A2 + A3 = and f5i + /32 + /^s = 0, which allow for some 
substitutions to be made. To avoid any ambiguity, we can rewrite them in terms of the fundamental weight 
basis ui = 5(2, —1, —1) and W2 = 5(1, 1, —2). Then if A = liL0i + l2t02 and /3 = bicoi+b2t02, the polynomials 
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take the form 







P3 


= n 




-2fei -62) 


P6 = 


n 


4(^1- 


^2^2 


-bi- 


262) 






P4 


= n 






Pi = 


n 




F2?2 


-61 4 








P5 


= n 


4(2^1- 


^^2 +&1 +262) 


P8 = 


n 


4(^1- 


^2^2 


+ 26i 


+ 62) 



The domains of polynomiahty of a weight diagram for a given A will be the (possibly empty) polytopes 
Ti nL(A) for i = 1, 2, . . . , 8, so there are at most eight of them, although in practice at most seven appear at 
a time. We could obtain a symbolic description of the domains of polynomiahty for the weight diagram of 
any A of A2 from the chamber complex. When A is one of the fundamental weights, we get a triangle with 
constant multiplicities inside; otherwise we get a hexagon with a (possibly empty) central triangle in which 
the multiplicities are constant and decrease linearly outside. Figure 0] shows what happens when we move 
from one fundamental weight to the other. The picture for A2 is already well-known (see, for example, 

[IHlEl). 




Figure 4: Weight diagrams and their domains of polynomiahty for A2. 



Corollary 12.61 explains why the second and third diagrams, as well as the fifth and sixth of Figure 0] are 
variations of each other, and why the polynomials attached to each of the seven regions are the same for 
each of these pairs of diagrams. For A2, the cones project under p\ to the three cones 



Co - pos((2,-l,-l),(l,l,-2)) = pos(a>i,c^2) 

Ci = pos((2, -1,-1), (1,0,-1)) = pos(wi,cji +W2) 

C2 = pos((l,0,-l),(l,l,-2)) = pos(wi+W2,^^2). 



We can see that Ci and C2 partition the fundamental Weyl chamber Co of A2, and hence c'f^ consists 
of Ci, C2 and all their faces. Therefore for A2 there are only two generic types of A's : A's with A2 < 
(diagrams 2 and 3 on Figure 0J and A's with A2 > (diagrams 5 and 6 on Figure EJ. The case A2 = 
corresponds to the regular hexagon, while the degenerate cases Ai = A2 and A2 = A3 correspond to the 
triangles. If we express A in terms of the fundamental weights A = litoi + ^2^2, these correspond to li < I2, 
h > I2 and li = I2 respectively for the hexagons, and Zi = and ^2 = for the degenerate cases (triangles). 
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We can write 
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7.2 The chamber complex for sUC {A3) 



with 



Ei = 



/ 





1 


1 


1 























\ 




1 


-1 


1 





1 




























1 











1 




























1 











1 






















-1 


-1 














1 
















-1 





-2 

















1 













-1 





-1 




















1 










-1 





-1 























1 





V 


1 


-1 





























1 / 



and B4 



P. 



( Ai + /3i + /32 + /33 
A2 — A3 
Ai — A2 
A2 — A3 
A2 + 2A3 - /3i - /32 - /33 
2A3 -fii-h 
A3-/3i 
A3-/32 
V A2 - /33 



Remark 7.1 E,^ is not unimodular. We don't know of a unimodular matrix for sl^C that would make the 
multiplicity function into a single partition function. 



For ^3 we must use the computer to do most of the computations. A symbohc calculator like Maple 
or Mathematica is especially useful. Here we used Maple (versions 7 and 8) and the package convex by 
Matthias Franz jSj. 

The set S^^^ of bases for i?4 has 146 elements, so there are 146 base cones for cr e B'-^\ These are 
9-dimensional cones, however they collapse to 132 6-dimensional cones when intersected with 

S = Q j : A e K*, /5 e M^^ Ai + ■ • • + A4 = 0, /3i + • ■ • + /34 = o| 



The full chamber complex is the complex of all intersections of the 146 cones in MP, and it has 6472 
full-dimensional cones. The chamber complex in i?-space has 1202 full-dimensional (6-dimensional) cones 
Tk, which we can rectify to get the chamber complex C^"*^ in (A, /3)-space with cones r^, k = 1, . . . , 1202. 
However, the 6-dimensional chamber complex thus obtained is not closed under the action of the symmetric 
group 64 on the /3-coordinates. 

Despite the fact that the chamber complex seems to lack the symmetry property in /3, we will see, as we 
find the polynomials attached to the domains of polynomiality, that there is a way to regain it. We can 
compute the polynomial associated to each of these 1202 cones by interpolation, for example using the fact 
that De Loera and Sturmfels computed the polynomials for the Kostant partition function for A3 in |25|. 
These 1202 polynomials are not all distinct. 

Observation 7.2 // we group together the top- dimensional cones from {rfe : k — I, . . . , 1202} with a 
particular polynomial, their union is always a convex polyhedral cone again. Grouping cones this way yields 
a glued chamber complex Q in (X, (3)-space with 612 cones Gk, k — \, . . . , 612. These cones form 64 orbits 
under the action of S4 on the P- coordinates. 



Proof. Here is a description of the algorithm used to make this observation. Suppose that {r^j , , . . . , Ti„ } 
consists of all the cones with a particular given associated polynomial, and let r be the convex polyhedral 
cone spanned by the union of all their rays. We want to prove r = uj^j^T^^. . 

We can find an affine half-space whose intersection with each of these cones is non-empty and bounded, 
so that we can work with truncated cones. The half-space Ai < 1 works. The union of {r^^, . . . ,Ti„} will 
equal r if and only if the union of their truncations gives the truncation of r. The truncated cones are 
polytopes, and we can compute their volume. We can check that the union of the truncations of r^^ , . . . , ri„ 
is the truncation of r just by checking that the volumes match. We know the Ti. have disjoint interiors 
because they are defined as the common refinement of base cones, hence the volume of the union of all 
these truncated cones is simply the sum of their volumes. If the computations are done symbolically (in 
Maple), there is no danger that tnmcated cones with very small volumes could create round-off errors. 
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The volumes are compared symbolically for every family of cones corresponding to the same polynomial. 
We glue together all the cones with the same polynomial, and observe that we still get a complex of convex 
polyhedral cones. This glued complex is invariant under the action of 54 on the /3-coordinates. I 



We now have two 6-dimensional chamber complexes C^^^ and Q, and we can construct the complexes Cj^ 
and Ga by first projecting all the cones through pA and then forming their common refinement. 

After transporting the hyperplane x + y + z + w = OofM* into the hyperplane z = through the orientation- 
preserving isometry 

/ 1 -1 -1 
T - 1 1-11 

V -1 -1 -1 

we can work in the coordinates (x, y, z) and look at the intersections of the complexes with the hyperplane 
z = l of R3. 

Figures |31[7| show the complexes intersected with the hyperplane z = 1 and also the complexes modulo the 
symmetry A ^-> — A''°^, which translates into a reflection along the central (vertical) line of the complexes. 
This symmetry reflects the symmetry of the Dynkin diagrams for An ■ Figure |H1 shows that even though 
regions appear and disappear along the lines of the complex (facets of the full-dimensional cones of the 
complexes), the complex given by simply looking at the number of regions in the permutahedra is coarser. 

Observation 7.3 For A^, only six generic cases occur. Generic permutahedra are always partitioned into 
213, 229, 261, 277, 325 or 337 regions. Degenerate cases occur along the walls in Figure\^ 



1\ 

-1 

-1 
-1 / 



(61) 



Projecting the cones of the glued complex Q on A-space gives 62 distinct cones, 60 of them corresponding 
to individual orbits under the action of ©4 on the /^-coordinates. The chamber complex Qa we get by 
taking their common refinement has 50 regions, or 25 modulo the symmetry A — A'''^^. This complex 
classifies the combinatorial types of A, i.e. the A's with the same partitioned permutahedra and family of 
polynomials. 




Figure 5: The chamber complex C 



We now give a simple proof of Theorem 13.21 for A3 and show that regions of the permutahedron given by 
(|23|l are as large as possible. 
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Figure 6: The chamber complex Qa_. 




[1.0,1] [0,0,1] 



Figure 7: The chamber complex Qa for Ai < — A4 in terms of fundamental weights. 



Theorem 7.4 For A3, the optimal partition of the permutahedron into domains of polynomiality for the 
weight multiplicities coincides with the partition of the permutahedron into domains of polynomiality for 
the Duistermaat-Heckman measure. 

Proof. We give a computer verified proof that for all A in the fundamental Weyl chamber of sUC the 
intersection of Q with L(A) defines walls within the permutahedron as determined by Theorem 13.11 We 
do this by expressing the walls of the permutahedron as the convex hulls of subsets of its vertices. The 
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following is an outline of our algorithm. 

The full dimensional cones Gi, . . . , Gei2 of the complex Q, when intersected L{X), subdivide the permutahe- 
dron into regions. For generic A, Gk ni(A) is either empty or a 3-dimensional region of the permutahedron. 
Furthermore, a 2-dimensional facet of that region will come from the intersection of a facet F of Gk with 
L(A), and an edge of that 2-dimensional facet will come from the intersection of a facet L oi F with L(A), 
and finally a vertex of that edge will come from the intersection of a facet of L with L(A). 

1. Set T equal to the set of all facets of the cones Gi, . . . , Gqi2- 

2. Classify the facets in JF according to their normals: call !Fi the subset of T consisting of all the facets 
with normal direction rii. Since each facet lies on a unique hyperplane, and since all these hyperplanes 
go through the origin, two facets will lie on the same hyperplane if and only if they have the same 
normals up to a scalar multiple. In our case, we find that there are 37 distinct normal directions. 

3. Set Ki ~ Upej^iF and verify that Ki is again a convex polyhedral cone. The verification is done 
by a truncation and volume comparison method similar to the one used in Observation 17.21 The 
intersection of the Ki with L(A) will be the walls partitioning the permutahedron. 

4. For each i, set Vi to be the set of facets of facets of Ki. The elements of V; are three dimensional 
cones. 

5. For each i, identify the / G Vi whose intersection with L(A) for generic A is a point. The convex hull 
of those points is Ki D L{X). These points are all vertices of the permutahedron, and the walls they 
define are exactly those of equation H23(l . 

We will illustrate this last step on an example. We find that one of the Vi consists of the 10 cones, which we 
will denote /i, . . . , /iq. One remarkable thing about the cones fj is that the first four coordinates of their 
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rays always correspond to one of the fundamental weights, while the last four correspond to a conjugate of 
the same fundamental weight. That is true for all V,. We have 



/l = 


pos(a;i,a;i), (a;2, tt ■ ^2), (^3,^3)) 


/6 


= pos((w2,7r- a;2), {w3,(f) ■ 0)3), {^3,(^3)) 


/2 = 


pos((a;i,wi), {w2,'JT-U2), {0)3,^-^3)) 


fr 


= pos((a;i, cr • wi), (wi,wi), (^3, • ^^3)) 


/3 = 


pos((a;i,c7-a;i), {u>2,tt ■ (J2), {u>s,(j) ■ 0J3)) 


fs 


= pos((a;i,a;i), {u3,(l) ■ 0)3), (^3,^3)) 


/4 = 


P0s{{u!i,a ■ iOi), {(J2,TT ■ (J2), (W3,W3) 


h 


= pos((a;i, cr • wi), (wi.wi), (ws.ws)) 


/5 = 


P0s{{wi,a-U)i), {U)3,4)-W3), {W3,Ul3)) 


fw 


= pos((a;i, cr • wi), {w2,n-W2)) 


where a = 


(1 3),7r=(2 3), </>=(2 4). 







To find the intersection of one of these cones with ^(A), we want to see whether there is a linear combination 
of its rays with nonnegative coefficients that would lie in L{X). If the rays are ri, . . . , rg, we are looking 
for ai , . . . , ttg > such that 



ai n + 02 r2 H h as = (Ai, A2, A3, A4, *, *, *, *) , 

or equivalently, 

aiPA(ri) + a2PA(r2) H \-PA{asrs) = (Ai, A2, A3, A4) , 

Hence we will get vertices for those A's and /j's such that A G Pa(/j )- So we compute the PA{fj)- 



Pa(/i) 


= pos(a;i,a;2,W3) 


PAifs) 


= pos(a;i,a;3) 




p0s(tJi, W2,t^3) 


PAifa) 


= pos (0)2,^^3) 




= pos(u;i, W2,t^3) 


PA{h) 


= p0s(u'i,W3) 


Pa{U) 


= pos(a;i,W2,W3) 


PA{h) 


= pos(ti;i, ws) 






PAifg) 


= pos(a;i,W3) 






Pa(/io) 


= pos(a;i,a;2) 



Only the first four of the cones span the fundamental Weyl chamber; the other six won't intersect L{X) for 
generic A. 

Observing that the last four coordinates of the rays of the /^ 's can always be obtained by applying a single 
permutation to the first four, we can rewrite /i, /2, /3, /4 as 

/i = pos((a;i,(2 3)-a;i), (a;2,(2 3)-a;2), (a;3,(2 3)-a;3)) 

/2 = pos((c^i,(2 4 3)-a;i), (0^2, (2 4 3) • c^a), (0^3, (2 4 3) ■ c^a)) 

fs = pos((u;i,(12 4 3)-wi), (w2, (1 2 4 3) ■ a.'2), (a.'3, (1 2 4 3) ■ ^3)) 

fi = pos((a;i,(12 3)-wi), (w2, (1 2 3) • W2), (^3, (1 2 3) • W3)) 



It then follows that 

/inL(A) = (A, (2 3) -A) 

/2nL(A) = (A, (2 4 3) -A) 

/3nL(A) = (A, (1243) -A) 

/4nL(A) = (A, (12 3) -A) 

which means there will be a wall with vertices 

(2 3) -A = (Ai,A3,A2,A4) = A' 

(12 3) -A = (A3,Ai,A2,A4) = (12)A' 

(243) -A = (Ai,A3,A4,A2) = (34)A' 

(1243) -A = (A3,Ai,A4,A2) = (12)(34)A' 



7.3 Further observations on A3 
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in the permutahedron for A. This waU is the convex huU of W ■ X' with W the parabohc subgroup 
((1 2), (3 4)). All these vertices have the same scalar product with (1, 1, —1, —1), and thus they lie on the 
same hyperplane with normal (1,1,— 1,-1). They are also the only points of ^4 • A lying on that hyperplane 
(for A generic). 

This process is automated, and we can express all the walls in this fashion by repeating this process for all 
the Vi's. We finally check that these walls are the same as those that partition the permutahedron for the 
DH-measure. ■ 

The set of domains for a given permutahedron is closed under the action of the Weyl group, so they come 
into orbits. Out of the 64 orbits of cones in the chamber complex Q, at least 22 orbits (when there are 213 
domains) and at most 31 orbits (when there are 337 domains) appear at a time. 



7.3 Further observations on 

With the chamber complex and all the weight polynomials for A3 in hand, we can test whether the bounds 
on the number of parallel linear factors from Theorems 16.21 and 16.51 are tight in this case. For k — 4, 
and a generic dominant weight, we should be getting at least two parallel linear factor in the directions 
conjugate to Wi or 0^3, and three parallel linear factors in the directions conjugate to lu2, for both the weight 
polynomials in the boundary regions and the jumps between adjacent regions. 

For generic A, two full dimensional domains are adjacent if the corresponding cones in the chamber complex 
are adjacent. This allows us to test for factorizations on the chamber complex level and thus work with all 
A at once, so to speak. We have verified that the bounds are met exactly here: when extra linear factors 
occur, they are not part of the parallel family 

For nongeneric A, we have to be careful when relating adjacency between regions of the permutahedron 
and adjacency between cones. This is because, in three dimensions for example, it is possible for two cones 
to touch along an edge but not along a face, and if we cut them with a hyperplane going through that 
edge, their two-dimensional intersections become adjacent because they have an edge in common. We get 
problems on the wall A2 = A3 of the fundamental Weyl chamber, for instance. The weight polynomials in 
this case all have the form 7(7 + 1), and the jumps between adjacent cones always have the form 

7(7 + 1) - 7' (7' + 1) - (7 - 7')(7 + 7' + 1) • 

Hence we don't get parallel linear factors at all in this case. For the cases Ai = A2 and A3 = A4, we get 49 
or 61 domains generically (see |14j for a study of the Duistermaat-Heckman function and its jumps for one 
of the 49 domain generic cases). The weight polynomials are at most quadratic in this case, and we verify 
that we get two parallel factors in every jump. 

The zero weight does not always appear in the weight space decomposition of A, but for A generic, there 
will be a non-empty domain of the permutahedron that contains the origin (even if it is not a weight) . This 
domain is invariant under the action of the Weyl group and we will call it the central domain. 

We can describe the generic central domains for A3 . The diagram on the left in Figure 1^1 shows the four 
types of domains for Ai < — A4. The light region corresponds to cubic central domains. In the region next 
to it, we get truncated cubes: four vertices in tetrahedra position in a cube are truncated, and we get six 
hexagonal faces and four triangular faces, forming a polytope with 16 vertices. In the remaining two darker 
regions, we get tetrahedra, in two different orientations (the central domain vanishes on the wall between 
them). When the hyperplane supporting a wall goes thought the origin, the bounded part giving the wall 
also contains the origin, and the central domain then vanishes (degenerates to a single point). This occurs 
when A = -A"^, A2 = or A3 = 0. 

The behavior of the weight polynomials, as polynomials in both A and f3, is slightly different. For generic 
A, we get only two of them (four if we don't work modulo the Dynkin diagram symmetry). The diagram on 
the right in Figure El shows the separation. The polynomial in the light region does not depend on the (3- 
coordinates, thus showing that permutahedra for A3 obtained by perturbing around the permutahedra for 
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the fundamental weights uji and lo^ have tetrahedral central domains over which the multiplicity function 
is constant (for fixed A). Domains like these are called lacunary domains in 14 . For Ai < — A4, the 
polynomials are 

- i(A2-A3 + l)(Ai-A2 + l)(Ai-A3 + 2) 

P^^"''^^ - ^(Ai - A2 + 1)(-A2 - 2X1 + A3A4 - A2A3 - A2A4 + A2 - A4 + 2 - 2/i2(/3i/32/33)) , 

where /12 is the complete homogeneous symmetric function: 

/»2(/3l, /?2, M - I3f + P2 + Pi + + 132133 + Plf33 . 




Figure 9: Central domains according to the shape of the polytope (left) and the weight polynomials (right). 

Remark 7.5 Dealing with A4 is more difficult computationally. For example, the permutahedron for the 
weight \ = 5 splits into 15230 regions, and this number is a lower bound on the number of maximal cells 
of the chamber complex for the weight multiplicities. 
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